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Abstract 



Criteria which a space-time must satisfy to represent a point mass embed- 
ded in an open Robertson-Walker (RW) universe are given. It is shown that 
McVittie's solution in the case k = satisfies these criteria, but does not 
in the case k = — 1. Existence of a solution for the case k = — 1 is proven 
and its representation in terms of an elliptic integral is given. The following 
properties of this and McVittie's k = solution are studied; uniqueness, be- 
haviour at future null infinity, recovery of the RW and Schwarzschild limits, 
compliance with energy conditions and the occurence of singularities. Exis- 
tence of solutions representing more general spherical objects embedded in a 
RW universe is also proven. 
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I. INTRODUCTION 



This paper deals with the embedding of massive objects in Robertson-Walker (RW) 
universes. There are three ways in which the physical embedding may be modelled and 
thus treated mathematically. Firstly, one can treat the body as a test body whose dynamics 
are described by a suitable set of equations of motion (e.g. geodesic equations for a test 
particle, equations derived from the Nambu action for cosmic strings). Secondly, the history 
of the surface of the object can be treated as a boundary S~, which is then matched with a 
diffeomorphic surface S + in the 'exterior' RW geometry. The usual matching conditions are 
continuity of the first and second fundamental forms of £ = S 1 * 1 |I| . This technique has been 
used to study the formation and evolution of voids in cosmology @, and to study domain 
walls §. 

The third method is to solve Einstein's field equations, exactly or approximately, in 
such a way that the resulting solution can be interpreted as an embedding of some massive 
object in a RW background. Two landmark papers in this vein are those of McVittie Q, who 
gave solutions of Einstein's equations with perfect fluid source which have been claimed to 
represent the embedding of the Schwarzschild field in the three families (k = 0, ±1) of RW 
space-times, and of Hawking ||, who studied gravitational radiation from a bound source 
in the k = — 1 dust filled RW space-time. There also exist several papers dealing with the 
superposition of the Kerr-Newman and RW space-times (see , or [0] for a summary) . We 
stress that we are considering the embedding of extended objects of finite size, so that the 
extensive studies (see e.g. ||) of perturbations which occur throughout the universe do not 
concern us here. 

These three approaches incorporate various degrees of coupling between the mass-energy 
of the extended body and the geometry of the universe at large. In particular physical 
situations, one of the three provides an appropriate model. For example, in cosmology, the 
galactic source of some observable effect is treated as a test body moving on a time-like 
geodesic of the RW geometry; the Einstein-Straus vacuole || provides a description of the 
effect of the cosmic expansion on the gravitational field of the sun, but McVittie's solution 
is a more appropriate description of the gravitational field outside a super-massive spherical 
body in an otherwise uniform RW space-time. 

An interesting result was derived recently by Senovilla and Vera regarding the cylindrical 
analogue of these models |10| . String dynamics in an RW universe deals with extended bodies 
which are limits of cylindrical objects and is well understood. However, moving to the next 
level of coupling, the aforementioned authors showed that no static cylindrical region can 
be matched continuously to a RW universe. Similiar results have also been obtained for the 
axially symmetric case |TT|. Since real strings have internal structure, this implies that at 
this level, strings cannot be embedded in a RW universe. The same is true for any static 
locally cylindrical objects; coins, bottles and (true) cylinders. This begs the question: Can 
the third type of embedding be carried out for cylindrical objects? i.e. can we find an 
exact solution of Einstein's equation representing a cylindrical object embedded in a RW 
space-time? 

We will not attempt to answer this question here, but by examining carefully the spherical 
case, suggest how the problem may be approached. Thus we readdress the problem first 
discussed by McVittie, but from a modern point of view. 
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We find that McVittie's solution in the case k = satisfactorily describes a massive 
particle embedded in a RW universe, but that his k = — 1 solution does not. Motivated by 
the differences between these two solutions, we lay down a priori conditions that a space- 
time (V, g) must satisfy to represent a massive particle embedded in a RW space-time. We 
provide a solution in the case k = — 1 and discuss uniqueness in each case, which has not been 
done before. We emphasise that we have not found a new solution of Einstein's equations, 
but have determined the physical significance of a certain class of shear-free spherically 
symmetric perfect fluid solutions (see chapter 14 of ||12|| ). Furthermore, we discuss the 
following properties (mathematical and physical) of the solutions: (a) representation of 
the k — — 1 solution by an elliptic integral, (b) recovery of the Schwarzschild solution in the 
vacuum limit, (c) behaviour at future null infinity, (d) compliance with energy conditions and 
(e) existence and nature of the central singularity. A central tool in this analysis is Hawking's 
quasi- local mass IJ. We show by this example how such quasi-local constructions can be 
used to obtain boundary conditions for Einstein's equations useful for obtaining solutions in 
particular situations. Since we are dealing with asymptotics in open space-times, the case 
k = +1 is excluded from our discussion. 

The structure of the paper is as follows. In the next section we review McVittie's solution 
and the Hawking mass and point out problems with the interpretation of the former in the 
case k = —1. In §111, we set out conditions for a space-time (V,g) to represent a massive 
particle embedded in a RW universe. Using these conditions and Einstein's field equations, 
we show how the problem reduces to finding a solution, with a certain asymptotic behaviour, 
of a non-linear second order differential equation. In §IV, we prove the existence of such 
a solution, and discuss uniqueness. Using this solution, the properties listed above are 
discussed in §V, and we make some concluding comments in §VI. The global struture of 
these space-times is to be analysed in a subsequent paper. 



II. MCVITTIE'S SOLUTION AND THE HAWKING MASS 

In 1933, McVittie || found solutions of Einstein's field equations for a perfect fluid 
energy-momentum tensor, representing a Schwarzschild field embedded in the RW space- 
times. His solutions can be written JL3] (using units in which c = G = 1) 




where 

m = m(t), (3 = P(t), (3 = -2(m/m) , (2.2) 

and here and throughout, a dot indicates partial differentiation with respect to t (a prime 
will be used for differentiation w.r.t. the variable r). The functions h(r),w(r) depend on a 
choice of k(= —1,0,-1-1), the Riemannian curvature of the surfaces of homogeneity t = const, 
in the background RW universe; 

{sinhr, k = — 1; f2sinh|, k = — 1; 

r, k = 0; w(r) = < r, k = 0; 
sinr, k = +1; [ 2sin~, k = +1. 
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The isotropic pressure p mv and the energy density p mv obtained from Einstein's field 
equations are given by 

q . .. / 1 _|_ JUL \ Up-P 
8^v = ^ 2 + I ^J F . (2.4) 



The properties of this solution have been summarized by Raychaudhuri as follows (cf. fli"4 



p. 97): "The McVittie solution follows uniquely under the following conditions: (i) The line 
element is spherically symmetric with a singularity at the centre, (ii) The energy-stress 
tensor is that of a perfect fluid, (iii) The fluid motion is shear free, (iv) The metric must 
asymptotically go over to the isotropic cosmological form." (It should be noted that neither 
a proof of this statement, nor a reference to one is given; McVittie's ad hoc approach does 
not include such a proof.) 

The function m = niQe' 13 ^ 2 for some constant m by (|2.2j ), and is interpreted as the mass 
at the singularity. When this is set equal to zero, the line element ( |2.1| ) is that of a RW 
space-time. 

The characterization of McVittie's solution quoted above is unsatisfactory, as points 
(i) and (iv) refer to properties which are deduced simply by looking at the metric tensor 
components relative to the line element ( J2.1[ ). We can show that point (iv) in particular 
is misleading. This point seems to imply that the solution corresponds to a point mass 
embedded in the RW geometry, so that the gravitational field is asymptotically that of a 
RW space-time. 

Hawking || has made this notion precise with a renormalized (against the RW back- 
ground) quasi-local mass measured at future null infinity X + . Since we are dealing with 
asymptotic regions of the space-time, we restrict our attention to the cases k = —1, 0. 



The Hawking mass is defined by analogy with the Bondi mass \T5\ of a bound source 
of gravitation in an asymptotically flat space-time; it measures the mass of a bound source 
of gravitation in an asymptotically RW universe. The additions to the total (infinite) mass 
from the RW background are subtracted away in a gauge invariant manner, as we describe 
now. The construction is valid in any space-time. 

We use the null tetrad {l a , n a , m a , rh a }, where l a is chosen to be an outgoing null vector, 
and take v to be an affine parameter along the integral curves of l a , so that 

Taking S to be a space-like 2-sphere orthogonal to l a and n a (so that v =constant on S), 
the quasi-local mass surrounded by S is defined to be 

M(S) = kJ (-$ 2 -a\ + $n + A) dS , (2.6) 

where 

, ■ 1/2 
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and the terms in the integrand have their usual meanings in Newman-Penrose notation. In 
the appropriate limits, M yields the Bondi mass and the ADM mass, and is the spherical 



version of Hayward's improved quasi-local mass JT6 



The renormalization is carried out by subtracting the local (fluid) matter which manifests 
itself in the Ricci tensor terms, and to leave the non-local gravitational terms. To do this, 
Hawking || has defined 

Mx(S) = K/($n + A) dS, (2.7) 

and 

M 2 (S) = k J(-y 2 -a\)dS . (2.8) 

In order that the bound source mass is measured at X + , 

M 2 = lim M 2 , (2.9) 

V—tOO 

is defined to be the mass of the model. 

Using a suitable null-tetrad, we can evaluate ( |2.8j ) for McVittie's space-time (notice that 
due to spherical symmetry A = o = 0). We find 

M 2 (S) = m Kr- (2.10) 

Thus for k — 0, wherein h = w = r, we have M 2 = m , and so the Hawking mass is 

M 2 = m , (2.11) 

which verifies the interpretation of mo as being the mass of a point particle embedded in the 
RW cosmos in this case. However in the case k — — 1, since r — > oo as v — > oo, we find that 

M 2 (S)^oo, as v^oo (2.12) 

and so the renormalized mass is infinite. Thus the k = — 1 McVittie solution does not 
represent a point mass embedded in a RW space-time. 

The main aim of this paper, then, is to provide a solution of Einstein's equation which 
does represent a point mass embedded in the k — — 1 RW space-time. 



III. DESCRIPTION OF SPACE-TIME REPRESENTING A POINT MASS 
EMBEDDED IN AN RW UNIVERSE 

In this section, we give three conditions (C1-C3) on space-time (V, g) which, if these 
conditions are satisfied, we postulate to represent a point mass embedded in a RW universe. 
These conditions are motivated by the discussion above. For convenience, we will use the 
symbol (M, g) to refer to such a space-time. 

Condition CI. (M,g) is spherically symmetric with a shear-free perfect fluid energy- 
momentum tensor. 
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Demanding a perfect fluid energy-momentum tensor and spherical symmetry are obvious 
requirements for the space-time we seek to describe; the requirement that the fluid flow lines 
be shear-free is a convenience that eases integration of the field equations without ruling out 
the existence of a solution. There does not seem to be any a priori reason why the shear 
should be set equal to zero. 



Under these assumptions, the line element can be written as 12 



ds 2 = -e u dt 2 + e» [dr 2 + h 2 {r)duj 2 } , (3.1) 

where v = u(r,t), // = fi(r,t) and h is an arbitrary function of r, which we may always 
assume is one of the functions h(r) of §11. Doing so maintains the connection with the 
corresponding forms of the RW space-times. 

The density and pressure obtained from Einstein's field equations are given by 

Svrp = -^ 2 e-» - e-'V' + + 2^// + 3^} , (3.2) 

8*J> = e~"{\w + v") + \v' 2 + \jW + //) + ^} - c -*{/i + \f - ^} . (3.3) 



The remaining field equations reduce to 

e " = fo-rt) , (3.4) 

VP - \^ - h -p! = F(r)e^ 2 , (3.5) 

where g(t) and F(r) are arbitrary functions of their arguments. 

The Weyl tensor is Petrov type D, and on a naturally occurring null tetrad, the only 
non-zero Newman-Penrose component is 

m 2 = - l F{r)e-^ 2 . (3.6) 
o 

Then the Hawking mass enclosed by any t = constant r = constant surface S is 

M 2 (S) = h 3 (r)F(r). (3.7) 
o 

This surface S is a metric sphere of radius R(r,t) := h(r)e^^ 2 . In the solutions we examine, 
dR/dr > 0, and so since M 2 (S) is independent of t, the limit r — > oo of (|3.7|) yields the 
Hawking mass at infinity. 

The second condition achieves two things. Firstly, it identifies the RW background 
against which the Hawking mass is measured, and secondly, it ensures that the limit r — > oo 
has physical significance. This condition is a minimal requirement that (M, g) 'looks like' a 
RW space-time near infinity. 

Condition C2. 

lim u(r,t) = (3(t) 

r— >oo 

for all t in the range of \i, and where (3 is the function appearing in the line element of the 
RW universe, 
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(3.8) 



As we shall see below, it is not necessary to give the corresponding condition for v [y — ► 
asr-> oo), as this is ensured by a naturally arising choice of function of integration. 

We turn now to the Hawking mass and consider how this may be used in our description. 
We wish to embed a finite non-zero mass in the RW universe, indicating that M 2 (S) should 
yield a finite positive number when measured at infinity, i.e. 

lim M 2 (S) = m (3.9) 

) — >oo 

for some positive constant rriQ. However this allows the possibility of surplus Hawking mass 
arbitrarily close to infinity. We will therefore consider the stronger condition, 

M 2 {S)=m , (3.10) 

for some positive constant itlq. This condition suggests that there is an isolated body of 
mass mo at the centre of the space-time. In order to present a more general discussion in 
§IV, we will deal with the condition (|3.9|), but focus on the special case of ( |3.1U| ). 
We see from the above that Q3.9|) implies 

F(r) = 0(h~ 3 ) as r — > oo , 

while (|3.10|) gives 



nr) = i^ y (3.1D 

This last equation is interesting, as it is a necessary and sufficient condition for the energy 
density to be spatially homogeneous, i.e. 

F(r) = constant x /T 3 (r) p = p(t) . (3.12) 

Thus to obtain the constant Hawking mass condition ( [3.11| ), or equivalently to express 
the fact that there is no 'extra' energy density in the universe outside the embedded mass, 
we take the third condition to be as follows. 

Condition C3. The energy density \3.2J obeys p = po(t), where po(t) is the energy 



density calculated via Einstein's field equations of the RW space-time with line element (3.8). 
Furthermore, the constant quantity m := h 3 (r)F(r)/6 is positive. 

We can show now how to obtain v —>■ oo as r — > oo. From the above, the energy-density 
of (M, g) is given by 

8np = Lrt) + 3ke -m 
4 

while that of the RW background is given by 
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(po and po, the energy- density and pressure of the RW background may be read off from 
( p.4|) and ( |2.3|) respectively by taking m = 0.) Thus p(t) = p (t) requires e 9 ^ = f3 2 . We 
will see below that there exists a solution satisfying the three conditions above for which p 
is differentiable with respect to t, and which obeys lim^oo p = (3. Thus 



lim e v = lim p 2 e 



-<?(*) 



as claimed. 

The physical problem has been modelled using the conditions C1-C3 above. We see 
that the remaining mathematical analysis is to find a solution p(r,t) of 

alii a 1 12 h /. 6m , , 

e M/ (H ~~k^> = ~hF> (3.13a) 



with the boundary condition 



lim //(r, t) = p(t) , (3.13b) 



for all t in the range of p. 

More generally, we look for a solution of (|3.13a|) with the right hand side replaced by 



0(h~ 3 (r)) as r — > oo. Our analysis below is based on this version of the equation, which 
corresponds to the condition (|3.9|). 

We note at this stage that McVittie's solution || in the case k — satisfies conditions 
C1-C3. We focus henceforth on the case k — — 1, and work under this assumption. 

We conclude this section by giving a useful transformation which will make the problem 
( |3.13| ) more manageable. Defining 7 = e~^l 2 and x = w 2 (c.f. §11), these become respectively 



7^ = G(x) 7 2 , (3.14a) 

lim 7 = e- /3/2 , (3.14b) 

where 

= -24m (x 2 + 4s) _5/2 = 0(x~ 5 ) as a; -> 00 , 

and the subscript indicates partial differentiation with respect to x. Using G(x) = 0(x~ 5 ) 
corresponds to the general case ( |3.9| ). We now proceed to prove existence of a solution of 
the boundary value problem ( |3.14|) . 



IV. EXISTENCE OF A SOLUTION AND UNIQUENESS CONSIDERATIONS 

Writing 7 (x, t) = a(t)(l + Y(x,t)) where a(t) := e~^l 2 allows us to restate the problem 
(OD as 

Y xx = aG + 2aGY + aGY 2 , 
Y(x, t) = o(l) , x — > 00 . 



S 



We treat x as a complex variable and solve the equation in a neighbourhood of infinity which 
includes {igC: Qx = 0, > x > 0} where x is some real constant. 
Consider the equation 

Y xx = aG + 5(2aGY + aGY 2 ), (4.1) 
where 5 > is a small parameter. We look for a solution of this equation of the form 

oo 

Y(x,t) = y £8 n Y n (x 1 t), (4.2) 

n=0 

and having found one, show that for sufficiently large values of x, this converges in the limit 
5 = 1 and obeys Y = o(l) as x — > oo. This will prove existence of the required solution. 
We will use the ansatz F(x, t) = a™ +1 (i)$/ n (x), and throughout this section, a prime on y n 
indicates differentiation with respect to argument. 



To proceed, we fill out ( |4.1| ) using (|4.2j), and equate powers of 5. This leads to the system 
of equations 



y>> = G(x) , (4.3a) 

n-1 

y" = 2G(x)y n - 1 + G(x) V m Vn~\- m , n>\. (4.3b) 



m=0 

Clearly, this system admits solutions obeying 

y n (x)=o(l), y' n (x) = o(x~ 1 ) , x -+ oo . 

Then we can write 

Vn(x) = f y' n (s)ds 



rx 

/ (x - s)y"(s) ds . 



Transforming the integral to one over a finite contour via z = s 1 and using basic bounds 
for integrals, we obtain 

\y n (x)\ < 2\x- x \ sup \s 3 y"(s)\ . 

{\x\<\s\} 

Then taking T> to be a neighbourhood of infinity contained in the intersection of {x G C : 
\x\ > x > 0} and a sector containing the positive real axis, we have 

\Vn{x)\ < C\\W n \\ V 

for all x E V, where C := 2xq , W n (x) := x 3 y"(x) and \\-\\x> is the supremum norm restricted 
to V. So 

\\yn\\v<C\W n \\v. (4.4) 



Applying these definitions to ( 4.3b|) , we obtain 
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n-1 

\W n \\ v < 2B\\y n ^\\ v + B II l/m||©||l/n-l-m||x' 

m=0 



where B = B(xq) = \\x G(x)\\-p, and so using (^4.4j ), 

n-1 

\\y n \\v < 2A\\y n _ 1 \\ v + A II 

m=0 

with A = BC. 

Using this inequality, we can derive a geometric bound for the ||yn|i>) which will suffice 
to prove the convergence properties of Vn required to show that this formal solution is a 
(convergent) solution. To see how, define the sequence of positive reals {& n }^o by 

bo = \\yo\\v , 

n—l 

bn 26 n _x -|- ^ * bra^n—l—m • 
m=0 

Then we see from the last inequality that 

\\Vn\v< A n b n ,n>0. (4.5) 
Consider the formal power series 

oo 

P(X) := J2 b n* n - 

n=0 

From the recurrence relation for the b n , we find that P(X) obeys 

P = b + X{2P + P 2 ) , 

the solution of which consistent with the definition of P is 

, , 1 - 2X - ((1 - 2X) 2 - AboX) 1 / 2 . s 

P{X) = " 2X ] ^_J_ ■ (4.6) 

This is an analytic function of X in a neighbourhood of the origin, and so for any < A G 1Z 
with |A| < radius of convergence of P(X), 

oo 
n=0 

is convergent. Then each term in this series must be bounded, i.e. there exists some positive 
real constant K such that 

b n \ n <K, n>l, 

and so 



W« £ K (a) 



A\ n 

n > 1. 
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Then for xq sufficiently large, this last inequality will read 



\\y n \v<K n , (4.7) 

for some < k < 1. To see this, notice that as Xq increases, the region D gets smaller, so that 
A = 2xq \x 3 G(x)\\x> decreases (recall that G(x) = 0(x~ 5 )). b = \\y n \\v is non-increasing, 
which by ( |4.6| ) indicates that the radius of convergence of P(X) in non-decreasing, allowing 
the use of non- decreasing values of A. 

The condition Q4.7Q is sufficient to imply that 

oo 
n=0 

converges uniformly on T> (see e.g. [0]). Hence 

oo 

Y(x,t) = Ea" +1 (%„(x) 

n=0 

converges uniformly on some subset of T> x TZ (x G T>, t G TZ) . By our construction, this is 
a solution of ( |4.1| ) in the case 5=1 obeying Y(x, t) = o(l) as x — > oo. 
Furthermore, for each fixed value of x, the series 

5>n(x)a" +1 (t), 5> n 0n)(n + l)a»(f)^ 

are both uniformly convergent on some interval of the real t-axis (which will contain the set 
{t G TZ : \a(t)\ < 1}). Hence by standard results dY/dt exists on this interval, and 

We summarise and extend as follows. 
Theorem 1 (Existence and uniqueness) There exists < Xq G TZ and a non-empty 
subset A C TZ such that on {x G TZ : |x| > Xo} x A, there exists a solution of r y xx = G{x) r ) 2 
where G(x) = 0(x~ 5 ) as x — > oo ; obeying j(x,t) = e~^' 2 + o(l) as x — > oo. T/iis solution 
is differentiate with respect to t. Furthermore, if G(x) is analytic in a neighbourhood of 
infinity, then this solution is the unique analytic solution. 

For the proof of the last statement, note first of all that if G(x) is analytic in a neigh- 
bourhood of infinity, then the existence proof of the theorem gives the construction of an 
analytic solution of ( p.l4j ). For each term in the series fl4.2p is found by integrating analytic 
functions, and is therefore analytic; uniform convergence guarantees analyticity of the sum. 
Next, we argue that if G(x) and a solution j(x,t) of (|3.14 ) are analytic functions of x in a 



neighbourhood of infinity, then this solution is unique. For let 71,72 be two such solutions, 
and define 

T := 7l - 72 , H(x, t) := ( 7l + l2 )G{x) = 0(x- 5 ) . 
Then H is analytic in x in a neighbourhood of infinity, and T obeys 
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Y xx = H(x,t)T, (4.8a) 
r = o(l),x->oo. (4.8b) 

By analyticity, the solutions of this linear differential equation, which has a regular singular 
point at infinity, can be written in the form 

oo 

r = J2 a n (t)x- n+p 

n=0 

for some real p. Using a Taylor series expansion about x~ l = for H(x, ■) and filling out 
the equation ( |4.8a| ), we find that the two independent solutions are described by 

(i) p = 0, ao arbitrary, a± = a 2 = 0, a n ,n > 3 determined by recurrence relations, and 
proportional to ao; 

(ii) p = 1, ao,ai arbitrary, a 2 = 0, a n ,n > 3 determined by recurrence relations given by 
linear combinations of ao and a\. 

In either case, we see that non-zero solutions do not obey V = o(l) as x — > oo, so that 
T = is the only solution of ([4.8|), and so 71 = 72, proving uniqueness. 



In particular, in the case of most interest to us, G(x) = — 24m (x 2 + 4x) _5//2 is such 
a function, and so the solution produced by the theorem is the unique analytic solution. 
Notice also that when itlq = in this solution, we obtain the RW line element ( |3.8j ). It will 
be useful to have the first few terms of this solution. These terms are obtained by integrating 
( |4.3a| ) with the appropriate choice of integration constants, and yield 

7 (x, t) = e~^{l - m( + 2 4x)1/2 - (x + 2) + (x 2 + Ax) 1 / 2 )} + 0{x- & ) , (4.9) 

where m = m§e~PI 2 . Converting to the original coordinates using the transformation given 
prior to ( |3.14j ), this leads to 

6 m = e /3(*)(i _ Am e- p/2 e~ 3r ) + 0(e~ 5r ) . (4.10) 
We note that the corresponding first order term in McVittie's k — — 1 solution is 

= e\\ + 2m e-^ 2 e- r / 2 ) + 0{e^/ 2 ) , 

from which we can identify the problem with this solution; the metric coefficients do not 
tend to those of the RW metric rapidly enough. 

We now turn our attention briefly to the case k = 0. The existence proof of this section 
is of course not needed for this case, because as we have seen already, McVittie's k = 
solution satisfies the conditions C1-C3. The uniqueness result of this section does apply, 
and so we see that if Raychauduri's conditions quoted in §11 are replaced by the conditions 
of §111, then under the added hypothesis of analyticity, the solution is indeed unique. 



V. PROPERTIES OF THE SOLUTION 

In this section, we discuss various properties of the unique analytic space-time (M, g) 
found in the previous section which obeys the conditions of §111 with k — — 1 and G(x) = 
— 24mo(a; 2 + 4a;) _5//2 . Throughout this section, terms such as 'the solution', 'the line element' 
etc. refer to this solution and its line element etc. unless otherwise specified. 
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A. Representation with an Elliptic Integral 



We show here that the function e M can be represented intrinsically by an elliptic integral; 
the basic results are from [|T2[] . 

Writing 7 = (x 2 + 4x) 1/,2 u, we can obtain the following first integral of ( p,14a| ); 



(x 2 + Axfu 2 x = A(u 2 - 4m u 3 ) + A(t) , (5.1) 

where A(t) is a function of integration. Using fl4.9|) above, we can compare powers of x -1 in 
this last equation to obtain 

A(t) = e-M . 
A straightforward calculation shows that 

(x +4xK = - T — , 

where R(r,t) := h(r)e^^ 2 is the radius of metric 2-spheres in the space-time, so to ensure 
dR/dr > 0, we take the negative square root of ( |5.1| ), which leads to an intrinsic elliptic 
integral representation for u; 

du 1 , / x \ „. . 

--In — - )+B{t), 



(4u 2 - 16m w 3 + e-Pyi 2 4 Vx + 4 

and -B(t) is the sole remaining function of integration. This term may involve mo, and so 
cannot be determined from the other relevant limit, m = The integral on the left hand 
side does not have a representation in terms of elementary functions for mo 7^ 0, and so nor 
does our solution. 

We note that it may be possible to determine the functions A(t),B(t) in the elliptic 
integral which yield the correct asymptotic behaviour for fi without prior knowledge of the 
solution. The advantage is that we would not need the existence proof of the previous 
section. However, this would be a rather difficult problem involving inversion of asymptotic 
formulae for elliptic integrals. We feel that the chosen method is the most direct, and has 
the advantage of dealing with the general case, G(x) = 0(x~ 5 ). 



B. Energy Conditions 



Using the field equation (|3.4j) and the first integral of the main equation (|3.13a|) as found 
in §V-A above, we can write the pressure ( j3.3|) as 



8vrj9 = + 3e- p - /T 1 /?^ + 2e~ p ) . (5.2) 

Notice that by the main theorem of §IV, in the limit r — > 00, this coincides with 87rpo 5 the 
pressure of the background RW universe. Thus all of the curvature tensor terms match up 
with those of the RW background in this limit. 
Using (fOD , we find that 



13 



8np = 8np + + 2 e - /3 )m e- /3/ V 3r + 0(e" 5r ) . (5.3) 
Recall also that the energy density obeys pit) = po(t). We see then that the question of 



whether or not p and p obey appropriate energy conditions |19| is, for sufficiently large values 
of r, equivalent to the same question regarding po and pq. Notice that if the weak energy 
condition is satisfied in the RW background, then p + p > 0, leading to (3 + 2e _/3 < 0. 
Thus according to (|5.3|) , to first order, the presence of a central mass in a k = — 1 RW 
universe causes a decrease in the fluid pressure, contrary to what one would expect. This is 
in distinction to the situation in the k = model, where according to (|2~3l) we can write 

8tt P = 8tc Po - (imoe-^r- 1 + 0{r~ 2 ) . 

The weak energy condition in the RW background implies (3 < 0, and so the first order 
perturbation of the pressure is positive, as expected. Indeed this behaviour is continued at 
all orders; 8irp > 87rpo i n the k = model, provided the weak energy condition holds in the 
RW background. 

This latter situation is in line with what happens in the analogous situation in Newtonian 
cosmology. The potential describing the physical scenario under consideration is found (in 
this linear theory) by adding the potentials m = — mr _1 and <p c = p(t)r 2 / 12 for respectively 
a point particle of mass m situated at r = and an isotropic cosmological model with density 
pit) (see [|20[ for the latter). We take the potential to be 

Tfl \ 

4> = 4>m + 4>c = H -77;p(t)r 2 ■ 

r 12 

Then the pressure across the surface {S : r = constant} at time t is 

P( s ) = / ~V0 • nd 2 S , 
J s 

where n is the unit inward normal to S. This yields 

2 

p(S) = Anm + -npr 3 . 

o 

Thus we see that in Newtonian theory, the central mass makes a positive contribution to 
the pressure. 

The negative first order contribution in the k = —1 case could be cancelled out by higher 
order terms, leaving a net positive contribution, but if not, it appears to be an interesting 
effect of the negative curvature of the spatial sections of the RW background. 



C. Recovering the Schwarzschild Space-Time 

We have seen above how the line element of the RW background is recovered by taking 
mo = 0. We show next how the line element of the exterior Schwarzschild field arises in a 
natural way as a limiting case of our solution. 

In McVittie's k = solution, the Schwarzschild field is obtained by setting $ equal to 
zero. Then following a constant rescaling of the coordinate r, the line element (|2.1| ) with 
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k = is the isotropic form of the Schwarzschild line element with mass parameter tuq. The 
procedure is quite natural; with (3 = 0, the energy density and pressure both vanish, yielding 
a spherical vacuum which is by necessity, the exterior Schwarzschild field. Note also that 
the expansion of the fluid flow lines 9 = |/3 is then also equal to zero. 
Carrying out the same procedure in the k = — 1 case leads to 

p(t)=p (t) = ^ 2 -3e^ = 0, 

giving 

= (t + c) 2 (5.4) 

for some constant c. Calculating Po(t), the pressure of the RW background, we find Po(t) = 
0, so that this space-time is (a portion of) Minkowski space-time. Similarly, calculating 
p(r, t) for our solution using this form of (3 yields p = 0, and so the Ricci tensor vanishes. 
From (|3.6|) , the Weyl tensor remains non-zero, and so by Birkhoff's theorem, (M,g) is (a 
portion of) the exterior Schwarzschild field. The Hawking mass of the Schwarzschild field 
is the Schwarzschild mass parameter, and so m in our solution is the Schwarzschild mass 
parameter. 

Thus McVittie's solution for k = 0, and our solution for the case k = — 1, represents the 
Schwarzschild field embedded in a RW universe. 

Another limiting case is of importance, namely when p + p = 0, so that the space-time 
is an Einstein space. In both cases (k = —1,0), we can explicitly verify that p + p = 
implies that p is constant. Again, the choice of f3(t) does not affect the value of the 
Weyl tensor, and so by the 'Birkhoff-with-a-cosmological-constant' theorem, space-time is 
a portion of the Schwarzschild-de Sitter cosmos. Thus the solutions discussed here give 
genuinely cosmological (i.e. non-stationary) generalizations of this static space-time. 



D. Behaviour at Future Null Infinity 

As r — > oo on the space-like hypersurfaces orthogonal to the fluid flow lines, the line 
element of our solution approaches that of a RW space-time. The question of how it behaves 
asymptotically along future null directions is more complicated, but the following argument 
indicates that the space-time tends to a RW universe in this limit. We show that the metric 
coefficients of our solution match up with those of the RW background as r — > oo along 
future null directions of the RW background, which are hence asymptotically future null 
directions of (M,g). We deal explicitly with the more complicated case k = — 1; analogous 
results hold for k = 0. We consider first the description of X + in the RW background. The 



following relies heavily on pi 



Consider the RW background, whose line element may be written 

rfs 2 = _ df 2 + e m | rfr 2 + sinh 2 rdu 2} 

= n 2 (r]) [-dr] 2 + dr 2 + sinh 2 rdu 2 } , (5.5) 

where VL{j]) = e /3 / 2 (t), dt = Q(r))dr] and duj 2 is the line element of the unit 2-sphere. Define 
coordinates u (0 < u < oo) and x (0 < X < °°) by 
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u = e v r <-> rj = — ln(u 2 + 2ux), (5.6) 
X = e v sinh r r = — ln(l + 2xu _1 ). (5-7) 



In these coordinates, the line element ( |5.5| ) assumes the form 

ds 2 = F 2 (u, x) { — du 2 — 2dudx + X 2 <^ 2 } ? 

where 

F(u, x) = ^ 2 (ln (u 2 + 2 MX ) 1 / 2 )( M 2 + 2UX)- 1 . 
Next, define Z = and introduce the non-physical line element 



ds 2 = H 2 (uJ)ds 2 

= -l 2 du 2 + 2M + du 2 , (5.8) 



where 



H(u, I) = 1ST 1 (hi (u 2 + 2ur 1 ) 1/2 )(u 2 + 2ur l ) l/2 . 



Then, in the usual way, future null infinity of the RW space-time is identified with the 
boundary H = of the space-time (V,g) whose metric is given via the line element ( |5.8|) . 
If H = coincides with I = 0, then a direct calculation shows that X + is a shear-free null 
hypersurface. This depends upon fi(^) being a sufficiently rapidly increasing function of 
its argument, which relates to the conditions required of e^(t) to ensure convergence of the 
solution in §IV. We note that for a perfect fluid with equation of state p = ap, we have 

/ 3a +1 \ 2 /( 3a+1 ) 
fl(rj) = A I sinh( — - — 77) 

for some constant A, which leads to 

H(u,l) = A _1 2 2//(3a+1) (5.9) 

l(u 2 l + 2u) . 

{{uH + 2m)( 3 «+ 1 )/ 2 - Z(3a+l)/2)2/(3«+l) ' 

so that these conditions are satisfied if 3a + 1 > 0. 

This shows how to describe X + in the RW backgrounds for a large class of such space- 
times. The importance of this for our situation is that it tells us that as x ~^ 00 along 
u =constant, we approach X + in the RW universe. To conclude this section, we simply 
note that in this limit, the metric coefficients of our solution approach those of the RW 
background, and all the curvature tensor terms approach the background values. Thus our 
solution is asymptotically RW at future null infinity; we already know it to be asymptotic 
to the RW background at spacelike infinity. 
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E. Singularities 



The solutions which we study here (McVit tie's k = solution and our k — — 1 version) 
have been shown to represent the Schwarzschild field embedded in a RW universe. It is 
therefore natural to ask if these solutions have a central singularity and event horizon, and 
if so, how they are affected by the cosmic expansion. We will treat this important issue in 
more depth elsewhere; we can give the following preliminary results here. 

For any space-time, the quantity / = ^o^ - 4^1^3+3^2 is an invariant of the curvature. 
Here, we have for both k = 0,-1 

Thus we see that there is indeed an intrinsic curvature singularity at the centre R = 0. The 
coordinates we have used might not cover this region; this is immaterial as (|5.11|) derives 
from an invariant property of the curvature tensor, namely M2(S) = itlq for all metric 
2-spheres S. 



Hayward |22| has shown that the Misner-Sharp gravitational energy is a useful tool for 
investigating singularities in spherical symmetry. One of the equivalent definitions for this 
quantity is 

E:=j(l- X ), X --=V a RV a R. (5.12) 



Carrying out a straightforward calculation which makes use of (|3.4j) and the first integral 
fl5.1|) , we obtain the following nice results, which apply to both k = and k = —1. 

X = - Y R2 p( f ) + 1 " SmoiT 1 , (5.13) 

4-7T 

E = — R 3 p(t)+m . (5.14) 

These forms have the advantage of being coordinate independent; both p and R are in- 
variantly defined quantities. ( |5.14|) is particularly satisfying; the effect on the gravitational 
energy of the presence of a particle of mass mo is an increase of exactly this amount. 
Notice now that if x is any point in the boundary R = 0, then 

lim E > , 

7— *x 

along any curve 7 approaching x. Thus by a result of Hayward P^| , the central singularity 
is space-like and trapped, as in the Schwarzschild space-time. 

In the case k = 0, we see from Q2.3p that there is also a curvature singularity at r = m/2, 
which, intriguingly, corresponds to R = 2mo, the gravitational radius of the central mass. 
We see from (|5.13|) that this is a space-like hypersurface, and is surrounded by a trapped 
region. The existence of this singularity is fundamentally different to the vacuum case, and 
demands a thorough investigation of the singularity and horizon structure of this space-time. 
These issues are currently being studied. 
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F. Summary 



The solution we have found represents a point mass m embedded in a k = — 1 RW 
universe. When m = 0, we obtain this RW background. The energy-density is identical 
to that of the background, and the zero-density limit gives Schwarzschild's space-time with 
mass parameter mo- The space-time is asymptotic to the RW universe at infinity and 
contains a space-like singularity at the centre. 



VI. COMMENTS 

We have given a prescription above for embedding the Schwarzschild field in a give open 
RW universe. Consider the converse problem. Given a spherically symmetric shear-free 
perfect fluid space-time {V,g), how do we know if (V,g) represents a point mass in a RW 
universe and if it does, how do we identify that RW universe? This presents us with a 
gauge problem. "Suppose we consider the lumpy universe model S, not knowing how the 
[background] model S was used to make the construction; can we uniquely recover S from 
S 1 ?" ||. In fact the answer to this question is yes. Calculate the Hawking mass for an 
arbitrary metric 2-sphere of (V,g). If the result is not a constant, then (V,g) does not 
represent a point mass in a RW universe. If the result is a constant (mo say), and if further 
when m = 0, (V,g) is a RW universe, we may proceed. This solves the gauge problem by 
identifying the background model. It remains then to check if (V, g) satisfies the remaining 
parts of conditions C2, C3 with respect to this well defined RW background. 

In a previous paper f23fl , we interpreted certain space-times as being extended sources 



for the McVittie field in the three cases k = 0, ±1. To further investigate the occurence of 
singularities and horizons in (M, g) , it would be interesting to determine if a collapsing fluid 
can be used as a source. This would allow us to interpret (M, g) as the end state of the 
spherical collapse of a massive body in an expanding universe and may throw some light 
on the issue of what, if any, the effect of this expansion is on the collapse. This still leaves 
the problem of whether the space-time is of black-hole (collapsed object surrounded by an 
event horizon) or white-hole (lagging core of an expanding universe) type. This issue is to be 
addressed in a subsequent paper in which the horizon, singularity and asymptotic structure 
of these space-times is analysed. 

The nature of the solution we have found has opened up these interesting questions. 
However, our main purpose was to give a clear physical interpretation of some solutions of 
Einstein's equation. In particular, we hope to have given such for McVittie's solution, which 
in the case k = does indeed represent a point mass in an RW universe; some authors have 



contested this interpretation We have seen how Hawking's mass was a useful tool 

in this. Our aim now is to use this tool in an attempt to identify solutions representing 
the embedding of other objects (cosmic strings, the Reissner-Nordstrom and Kerr fields) in 
RW universes. We note that some of the solutions given previously which it was claimed 
represent such do not reproduce McVittie's solution in the k = case or our solution in 
the k = —1 case in the appropriate limit (charge-free and non-rotating). This may be an 
inherent discontinuous feature of solutions of the field equations in such situations. However 
it leads to the suspicion that these solutions do not satisfy a set of conditions analogous to 
CI — C3 which clearly determine their physical interpretation. 
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